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MASSACHUSETTS 


An  Analytic  theory  for  the  evolution  of  high  altitude  satellite  orbits  is 
developed  in  this  note.  The  distinctive  feature  of  the  theory  lies  in  the 
double  averaging  of  the  differential  equations — once  over  the  period  of  the 
orbit,  and  secondly  over  the  period  of  the  moon.  This  technique  is  called 
"intermediate  averaging"  to  distinguish  it  frcm  the  conventional  doubly 
averaged  theories,  and  to  denote  the  time  scales  inherent  in  the  averaging 
technique. 
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I . INTRODUCTION 


An  analytic  theory  for  the  evolution  of  high  altitude  satellite  orbits  is 
developed  in  this  note.  The  distinctive  feature  of  the  theory  lies  in  the 
double  averaging  of  the  differential  equations — once  over  the  period  of  the 
orbit,  and  secondly  over  the  period  of  the  moon.  This  technique  is  called 
"intermediate  averaging"  to  distinguish  it  from  the  conventional  doubly 
averaged  theories,  and  to  denote  the  time  scales  inherent  in  the  averaging 
technique. 

II.  BACKGROUND 

Analytical  theories  for  the  evolution  of  orbits  of  artificial  satellites 
of  the  earth  have  been  developed  from  the  early  days  of  satellites.  The  aim 
has  been  to  provide  computationally  efficient  formulations  that  could  be  used 
' • bdtn  fbt*  orbital ‘analysis- and  T>rbit* estimation.  - (A  -comprehensive- bibliography-., 
is  provided  by  Dallas:  Ref.  1).  The  alternate  technique  is  the  accurate, 
relatively  compact  but  computationally  expensive  method  of  direct  numerical 
integration  of  the  differential  equations  for  orbital  evolution. 

By  far  the  most  common  technique  that  has  been  applied  in  the  develop- 
ment of  analytical  theories  is  that  of  averaging  over  the  time  scales  in- 
herent in  the  differential  equations.  An  adequate  exposition  of  the  method 
of  averages  is  provided  by  Lorell  et  al.  (Ref.  2) . The  characteristic  features 
of  orbital  theories  developed  using  the  technique  are: 

1.  They  are  expressed  as  analytical  formulae  for  the  rates  of  change  of 
mean  elements. 
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expansion  of  the  geopotential  (higher  order  theories  do  exist  and  are 


3.  The  first  averaging  is  over  the  mean  anomaly  of  the  satellite  which 


is  the  fastest  periodic  variable  in  the  system,  with  the  necessary 


short-periodic  terms  also  being  formulated 


which  is  typically  the  second  fastest  variable,  at  least  for  low 


A complete  theory  would  then  consist  of  the  doubly  averaged  secular  rates  of 


change  of  the  mean  elements,  along  with  the  long-periodic  and  the  short 


periQdip  tqppis 


III.  THE  SYSTEM  OF  EQUATIONS 


The  well-known  Gauss  form  of  the  variational  equations  in  the  Keplerian 


elements  of  the  orbit  are  used  in  this  paper  to  develop  an  alternate  theory 


for  high  altitude  satellite  orbits.  The  equations  are  as  follows 


sin  v + F 


and 


-jt7  = - (f1)  [-F  cos  v + F sin  v(l  + -)  ] - (cos  I)  — 

dt  e v 1 2 p dt 


p = ae 

V = the  gravitational  constant  of  the  earth 
r = the  geocentric  radial  distance  to  the  satellite 
v = the  true  anomaly  of  the  satellite 
n = the  mean  motion 

F , Fj,  F3  = the  components  of  the  perturbing  force  along  the  instan- 
taneous radial,  transverse  and  normal  directions  to  the 
satellite  orbit. 

Certain  features  of  the  above  equations  should  be  noticed.  Firstly,  the 
Gauss  form  has  been  altered  slightly  in  the  equations  for  e,  u and  M in  order 
to  take  advantage  of  commonality  of  terms.  Secondly,  the  differential  equa- 
tions are  linear  in  the  perturbing  forces.  Thirdly,  singularities  at  e = 0 
or  I ■ 0 exist;  these  will  be  commented  on  later. 


The  perturbations  considered  here  result  from  the  zonal  harmonic  terms  in 
the  geopotential  and  the  point  mass  gravitational  effects  of  the  sun  and  the 
moon.  The  perturbing  force  due  to  the  zonal  harmonics  can  be  written  as 


n = the  degree  of  the  zonal  harmonic 


the  zonal  harmonic  coefficient 


the  n-th  order  Legendre  polynomial  in  the  argument 


d(sin  6)  n 


thV “ifis\anl6otiS ' unit* vectors -a long-  the  radial,,  tj^asyerse.  and. 


normal  directions  in  orbit 


The  perturbing  force  due  to  a third  body  (the  moon  or  the  sun)  can  be 


the  order  of  the  term  in  the  Lengendre  polynomial  expansion  of  the 


the  q-th  order  Lengendre  polynomial  in  the  argument  z,  = (u  . u.  ) 


the  instantaneous  geocentric  unit  vector  to  the  third  body 


Geopotential  zonal  harmonic  terms  up  to  order  4 only  will  be  included  in 
the  theory.  Only  the  first  term  in  the  expansion  of  the  third  body  perturba- 
tion will  be  considered.  Neglected  terms  are  at  least  an  order  of  magnitude 
smaller,  for  satellites  in  half-synchronous  orbits,  and  at  least  a factor  of 
4 smaller  for  synchronous  satellites.  Atmospheric  drag  is  neglected  in  this 
theory. 

V.  THE  SINGLY  AVERAGED  EQUATIONS 

The  first  step  in  developing  an  analytical  theory  is  to  average  the  dif- 
ferential equations  over  the  fastest  variable  on  the  right  hand  side  of  the 
equations,  which  is  the  mean  anomaly  of  the  satellite.  A considerable  body 
of  results  already  exists  for  this  step  and  we  will  adopt  the  results  of  Liu 
(Ref.  3)  for  the  singly  averaged  equations  due  to  the  geopotential.  Liu 
tabulates  both  the  averaged  rates  and  the  short  periodic  terns*,  the  latter* 
being  computed  only  to  order  J^. 

The  third-body  perturbations  must  also  be  averaged  over  the  period  of 
the  satellite.  The  relevant  results  for  the  averaged  effects  are  taken  from 
earlier  work  (Ref.  4)  of  an  author  of  this  note.  The  short-periodic  terms 
due  to  the  third  body  were,  however,  not  derived  in  Ref.  4,  and  hence  were 
developed  for  this  note.  They  are  tabulated  in  Appendix  1 along  with  the 
singly  averaged  rates  of  the  elements  due  to  the  third  body. 

VI . INTERMEDIATE  AVERAGING 

The  singly  averaged  equations  due  to  the  geopotential  show  a periodic 
dependence  on  the  argument  of  perigee.  However,  the  rate  of  change  of  the 
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argument  of  perigee  in  high  altitude  orbits  is  less  than  1 deg/day.  Thus, 
the  period  is  of  the  order  of  a year  or  larger. 

The  singly  averaged  equations  due  to  the  third  body  exhibit  dependence 
on  the  position  of  the  third  body,  and  hence  on  its  period.  And,  for  the 
moon,  the  period  is  29  days.  Hence,  the  fastest  variable  left  in  the  singly 
averaged  equations  is  the  mean  anomaly  of  the  moon.  Mathematically, 

"If"  = fG(*'  jn)  + V*  V 

where  < > = singly  averaged  rates 

f_  = terms  due  to  the  geopotential 
G 

f = terms  due  to  the  third  body 

Mt  = the  mean  anomaly  of  the  third  body 

z = any  element  of  the  satellite  orbit  in  the  set  z. 

“ ' * • • - ....... 

The  next  step  in  averaging  is 

"at”  ■ V-  V * 7-  £"  V"  Vdt 

M */0 

where  t is  the  moon's  period.  Neither  the  geopotential  terms  nor  the  third 
M 

body  terms  due  to  the  sun  are  affected  by  this  step  of  averaging.  But  the 
mean  anomaly  of  the  moon  is  averaged  out  of  the  equations.  The  rates  aver- 
aged over  the  intermediate  period  for  the  moon  are  tabulated  in  Appendix  2. 

The  theory  would  not  be  complete  without  the  intermediate  periodic 
terms,  whose  basic  period  is  that  of  the  moon.  These  terms  have  been  derived 
and  are  also  given  in  Appendix  2. 

VII.  PROBLEM  OF  SINGULARITIES 

The  original  differential  equations  are  singular  at  the  values  e = 0 
and  1=0.  These  singularities  do  carry  over  into  the  theory  developed  here. 
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As  the  eccentricity  goes  to  zero,  the  perigee  becomes  ill-defined.  Hence, 
the  short-periodic  corrections  in  the  mean  anomaly  and  in  the  argument  of 
perigee  become  large.  However,  these  two  corrections  are  of  nearly  equal  mag- 
nitude and  opposite  in  sign.  Thus,  the  position  of  the  satellite  in  the  orbit 
as  defined  by  id  + M does  not  exhibit  any  singularity.  Unless  the  eccentricity 
becomes  so  small  that  the  precision  of  computation  affects  results  (an  un- 
likely case),  the  theory  remains  valid. 

As  the  inclination  goes  to  zero,  the  line  of  nodes,  and  hence  the  value 
of  ft,  become  ill-defined.  The  theory  does  exhibit  this  singularity  in  the 
ftcular  terms.  And  again,  the  change  in  the  value  of  ft  is  nearly  equal  and 
opposite  in  sign  to  that  of  id,  the  argument  of  perigee.  Thus  the  position 
of  perigee  can  be  defined  using  the  value  of  (ft  + id) . Another  solution  is 
also  possible:  the  reference  plane  can  be  changed  from  the  equator  to,  say, 
the  ecliptic,  or  preferably,  a plane  containing  the  earth's  polar  axis.  Very 
little  needs  to  be  changed  in  the  theory  to  accommodate  a plane  change. 

Finally,  the  classic  problem  of  singularity  at  critical  inclination  does 
not  exist  in  this  theory  as  no  averaging  is  carried  out  over  the  period  of  the 
argument  of  perigee. 

VIII.  CONCLUSION 

A theory  for  high  altitude  orbits  has  been  developed  using  the  concept 
of  intermediate  averaging.  The  final  theory  consists  of  the  following  addi- 
tive, separable  parts: 

1)  The  averaged  rates  of  change  of  the  orbital  elements  due  to  the  geo- 
potential (intermediate  averaging  leaves  the  singly  averaged  expres- 
sions unchanged) ; 
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2)  The  short  periodics  due  to  the  geopotential; 

3)  The  intermediate  averaged  rates  due  to  the  moon; 


* 


* 
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4)  The  intermediate  periodic  terms  due  to  the  moon; 

5)  The  averaged  rates  due  to  the  sun  (intermediate  averaging  leaves  the 
singly  averaged  expressions  unchanged) ; 

• I 

6)  The  short  periodic  terms  due  to  the  sun  and  the  moon. 

The  theory  is  being  implemented  in  an  analytic  orbit  determination  program 

« { 

called  ANODE  (Ref.  5)  for  the  Millstone  Hill  Radar. 


| | 
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APPENDIX  1 


SINGLY  AVERAGED  RATES  AND  SHORT  PERIODICS  DUE  TO  THIRD  BODY 

The  singly  averaged  rates  of  change  of  the  orbital  elements  of  a satellite 
due  to  the  point  mass  perturbations  of  a third  body  are  recorded  in  this  part 
of  the  appendix.  The  theory  is  restricted  to  the  first  term  in  the  Lengendre 
polynomial  expansion  of  the  perturbing  force  due  to  the  third  body, 
da 

— Q 

dt 

<||'>  = f l <5  ~ 4e,^i^3  cos  w “ e^2^3  Sin 


<d£>  = 3 A_  1 

dt  2 L sin  I 
e 1 


[(5  - 4eK1E;3  sin  u>  + e?2^3  cos 


<— > = Aee^E  E 

dt  2 ^1^2 


^ = I Ae>5(4^i  ~ ^ ‘ l)  " COS  1 <f^> 


cdM 

dt 


n>  = £ [8  + 12e  + 15e2  - 3(1  + 12e  + 22<*2)^  - 21e£2] 


- e^cos  I <§>  + <^>) 


where 


. , k,  .a  . 3 

A = n(— ) (-) 

k 

n = the  mean  motion  of  satellite 


y = gravitational  constants  of  k-th  body  and  the  earth 
a = semi-major  axis  of  satellite  orbit 
r^  = geocentric  distance  to  third  body 
I = inclination  of  the  satellite  orbit  to  equator 


e = eccentricity  of  satellite  orbit 

, 2 
e = 1 - e 

to  * argument  of  perigee  of  satellite  orbit 
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5 , £ , f = projections  of  the  geocentric  unit  vector  to  the  third  body 
12  3 

— ► ' -v 

on  the  satellite  orbital  axes  u , unit  normal,  u , to 

n P 

perigee,  and  u^,  normal  to  up  in  the  orbit  plane 
M = mean  anomaly  of  satellite. 

The  short-periodic  perturbations  due  to  a third  body  are  tabulated  below. 


A I = 3 - — t-  [£,£,  cos  oi(— ~ e sin  E + ■1~2E-  + 7—  sin  2E  - eJ  sin  E 
p n ^*3  1 3 4 4 2 


e 5e  1 e 

- — sin  3E)  + e sin  w(-^—  cos  E + — cos  2E  + — cos  2E 

■ 1 i*  5e  ^ g 

- — cos  3E)  + e ^2^3  cos  cos  E “ ^ cos  2E  _ 4“  cos  2E 

+ cos  3E)  - e sin  u)(-§-  sin  E - 7 sin  2E  + — sin  3E)  ] 


Aft  = 3 - — j— ^ [eVc,.  cos  (d(|^  cos  E - 7 cos  2E  - f-  cos  2E 

pn*j._ij  4 4 4 

c e sm  1 

+ 77  cos  3E)  + sin  w(-~  e sin  E - e3 *  sin  E + 7 sin  2E 

12  134  4 

fp-  e 1,  5e  1 

+ — sin  2E  - fr  sin  3E)  + e LL  sin  u(—  cos  E - 7 cos  2E 


^ sl 

7-  cos  2E  + vv  cos  3E)  + eC-C,  cos  w(-j  sin  E --  sin  2E 


+ Yj  sin  3E)  1 

Ae  - (Aa  ) - Be*5  7 7 [e*<d  - d)  (f  e cos  E - j cos  2E 
p 2ac  p n e 2 14  4 

2 1 e2 

- cos  2E  + cos  3E)  + ^1^2  (I  Sin  2E  ~ 2e  sin  E " sin  2E) 

2 2 

- eC, £„  (4  sin  E + ^7-  sin  E + 7-  sin  3E  - 77  sin  3E  - 7 sin  2E) ] 

12  2 4 b 2 


Au>p  = - cos  I(Aftp)  + ^ “ [«  sin  E - j sin  2E)  (1  - 3£2) 

- £ (-cos  E + | cos  2E)J  + 3-  - [£,£.,  (7-  sin  2E  - 2e  sin  E) 

i *■  4 ne!24 


+ Y~  (?2  “£]_)  sin  E + ~2  sin  2E  ~ g sin  3E*  + — 4 — e s*n  2E1 
AM  = -e^lcos  I(AQ  ) + (Am  ) ] - 2— [3e(sin  E + 7—  sin  E - 7-  sin  2E 


+ sin  3E)  + 3£2(-  ^ e sin  E - e3  sin  E + ^ sin  2E 


+ y~  sin  2E  - YJ  sin  3E)  + ^ ^2e^_  f sin  E “ J sin  2E 

e ij  5s  1 

+ T sin  3E)  + 6£  £ e (—  cos  E - 7-  cos  2E  - 7—  cos  2E 
6 1 2 4 4 4 


+ Y2  cos  3E) 


Aap  = ^ ^ [ec^OC^  ~ D (-cos  E + j cos  2E)  + 3eZ^Z2(S  sin  E 

- | sin  2E)  + 3^a22  - £2)  (e  cos  E - COS42E)  + 3C1?2  (J  sin  2E 


2e  sin  E)  + — e£  £2  sin  2E] 


where  E is  the  eccentric  anomaly  of  the  satellite,  and  6 = 1 + e . 


APPENDIX  2 


INTERMEDIATE  AVERAGING  OF  THE  PERTURBATIONS  DUE  TO  THE  MOON 


This  part  of  the  appendix  records  the  intermediate  averaged  perturbations 
of  a satellite  orbit  due  to  the  point  mass  effect  of  the  moon. 

<<4—>>  = -t  — - B sin  J,  cos  J,  1(5  - 4e)  cos  oo  sin  AP, 
dt  4 k k k 

e 

- e sin  a)  cos  AP^] 

dft  3 n B _ _ tic  . . . 

«— » = ; : — sin  J.  cos  J.  [(5  - 4e)  sin  u>  sin  AP, 

dt  4 *5  sin  I k k k 

E 

+ e cos  a)  cos  APjJ 

<<4r>:>  = nBee*3  sin  2AP,  sin2  J 
at  o K K 

«4r->>  = 4 nBe^lsin2  J,  cos  2AP  - 2]  - cos  I <<§4>> 
at  4 K k at 

- n>>  = -e^lcos  I «4r»  + «Gr»]  + x B[2(8  + 12e  + 15e2) 

at  at  at  o 


2 2 2 
- 21e(sin  AP^  + cos  cos  AP^) 


2 2 2 2 
- 3(1  + 12e  + 22e  ) (cos  AP^  + cos  sin  AP^) ] 


«n»  - » 

Intermediate  averaging  also  yields  periodic  terms  with  the  moon's  period. 

These  are  recorded  below. 

Aa  =0 
I,P 


* ■ . • 


AI  p = jjJ-^r{(5-  4e>  cos  0)  x (T)-  e Sin  aj  x (5)  } 
nJc  e 

® - (j  cos  AP^  cos  (2^  + 2v^)  + ^ sin  cos  sin  AP^  sin 

ek  6k 

+ — cos  Apk  cos  (2\  + 3v  ) + — cos  AP^  cos  (2(1)^  + v^) 


ek 

>k  cos  (2u^  + 3v^)  + — cos  AP^  cos  (2(1^  + v^) 

ek 

\ sin  J,  cos  J.  sin  v,  + — sin  AP,  sin  J,  cos 
k k k k 12  k k 


(2oj^  + 2v 


_ % 


k k 

■ — sin  AP^  sin  cos  sin  + — sin  APk  sin  Jk  cos  Jk 

ek 

iin  (20)^  + 3\>^)  + — sin  AP^  sin  cos  sin  (2^  + v^)  } 

sin  AP^  cos  (2(1)^  + 2v^)  + ^ cos  AP^  sin  cos  sin  (2^ 

e,  e, 

)c  )c 

— sin  AP^  cos  (2^  + 3v^)  - — sin  AP^  cos  (2oi^  + v^) 


+ 2V 


ek 

cos  AP^  sin  Jj 


cos  si 


ek 

Fk  sin  \ + 12 


ek 

Yj  cos  APk  sin  cos 


r,  cos  j, 
k ] 


ek 

(2<i>k  + 3\>k)  + — cos  AP^  sin  cos 

= [ (5  - 4e)  sin  w x iT}  + e 

2 nk 


fk  sin  (2wk  + V 

>s  a)  x 12,  ] 


are  as 


above  in  AI^ 


6 

- r Sr  Bechli sin  2AFk  {_ek sin  vk  - t sin  (2\  + V 

k 

1 ek  1 2 

- — sin  2(o>k  + vk)  - — sin  (2u>k  + 3v^)  } + — sin  2AP^  cos 

ek  i 

{ek  sin  vk  - j~  sin  (2«i>k  + vk)  - ~ sin  (2a»k  + 2vk) 

ek  1 ek 

- sin  (2o),  + 3vu) } + f cos  2AP,_  cos  {-  — cos  (2o>^  + v,.) 


where 


Aw 


I,P 


- j cos  (2wk 

3 n a,.1*,., 

- — Be  [4  x 
2 \ 


+ 2\>k)  - ^ cos  (2u>k  + 3vk>}] 

3 - , 4 ' - 5 ) - cos  I Aft 

' I , P 


^3,  = {j  cos2  APk  sin  (2a>k  + 2vk)  - ^ sin2  APk  cos2  Jk 

sin  (2wk  + 2vk>  - ^ sin  2APk  cos  Jk  cos  (2wk  + 2\>k) 

ek  2 ek  2 . 

+ — cos  APk  sin  vk  + — cos  APk  sin  (2wk  + 3vk> 

ek  2 ek  2 2 

+ — cos  APk  sin  (2wk  + vk>  + — sin  APk  cos  Jk  sin  vk 

Sk  2 2 

- — sin  APk  cos  Jk  s^n  ^2a,k  + 3\^ 

ek  2 2 ek 

- — sin  APk  cos  Jk  sin  (2wk  + \>k)  - — sin  2APk 

cos  Jk  cos  (2wk  + 3vk>  - g—  sin  2APk  cos  Jk  cos  (2wk  + vk> } 

12  12  2 
4)  = {^-  sin  APk  sin  (2wk  + 2\>k)  - — cos  APk  cos  Jk 

sin  (2^  + 2vk>  + ^ sin  2APk  cos  Jk  cos  (2wk  + 2vk> 

6k  2 6k  2 

+ — sin  APk  sin  vk  + sin  APk  sin  (2wk  + 3vk> 

ek  2 6k  2 2 

+ — sin  APk  sin  (2u>k  + vk)  + — cos  APk  cos  Jk  sin 

ek  2 2 ek  2 2 

- — cos  AP.  cos  J.  sin  (2w.  + 3\>.)  - j-  cos  AP.  cos  J 


sin  (2u>  + v ) + 


e,  sin  v, 


mean  motion  of  moon  around  earth 


parameter  of  the  orbit  of  the  moon 


inclination  of  the  satellite  orbit  with  respect  to  the  moon’s  orbit 
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